String Cosmology in Anisotropic Bianchi-II Space-time 



Suresh Kumar 



Department of Applied Mathematics, 
Delhi Technological University (Formerly Delhi College of Engineering) , 
Bawana Road, Delhi-110 042, India. E-mail: sukuyd@gmail.com 

Abstract 

The present study deals with a spatially homogeneous and anisotropic Bianchi-II cosmological model 
representing massive strings. The energy-momentum tensor, as formulated by Letelier (1983), has been used 
to construct a massive string cosmological model for which the expansion scalar is proportional to one of 
the components of shear tensor. The Einstein's field equations have been solved by applying a variation 
law for generalized Hubble's parameter that yields a constant value of deceleration parameter in Bianchi-II 
space-time. A comparative study of accelerating and decelerating modes of the evolution of universe has 
been carried out in the presence of string scenario. The study reveals that massive strings dominate the 
early Universe. The strings eventually disappear from the Universe for sufficiently large times, which is in 
agreement with the current astronomical observations. 
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£jy 1 Introduction 

In recent years, there has been considerable interest in string cosmology. Cosmic strings are topologically stable 
objects which might be found during a phase transition in the early universe (Kibble [1]). Cosmic strings play 
an important role in the study of the early Universe. These arise during the phase transition after the big bang 
qq ' explosion as the temperature goes down below some critical temperature as predicted by grand unified theories 
(Zel'dovich et al. [2]; Kibble [313]; Everett g]; Vilenkin [S]). It is believed that cosmic strings give rise to density 
perturbations which lead to the formation of galaxies (Zel'dovich [5]). However, recent observations suggest that 
cosmic strings cannot be wholly responsible for either the CMB fluctuations or the observed clustering of galaxies 
t-H ' The cosmic strings have stress-energy, and couple to the gravitational field. Therefore, it is interesting 

to study the gravitational effects that arise from strings. The pioneering work in the formulation of the energy- 
momentum tensor for classical massive strings was done by Letelier [5] who considered the massive strings to 
be formed by geometric strings with particle attached along its extension. Letelier |10j first used this idea in 
obtaining cosmological solutions in Bianchi-I and Kantowski-Sachs space-times. Stachel has studied massive 
string. 

The present day Universe is satisfactorily described by homogeneous and isotropic models given by the 
FRW space-time. But at smaller scales, the Universe is neither homogeneous and isotropic nor do we expect the 
Universe in its early stages to have these properties. Homogeneous and anisotropic cosmological models have been 
widely studied in the framework of general relativity in the search of a realistic picture of the Universe in its early 
stages. Although these are more restricted than the inhomogeneous models which explain a number of observed 
phenomena quite satisfactorily. A spatially homogeneous Bianchi model necessarily has a three-dimensional 
group, which acts simply transitively on space-like three-dimensional orbits. Here we confine ourselves to models 
of Bianchi-II. Asseo and Sol [T2] emphasized the importance of Bianchi type-II Universe. Bianchi type-II space- 
time has a fundamental role in constructing cosmological models suitable for describing the early stages of 
evolution of Universe. 

Roy and Banerjee [13] have dealt with locally rotationally symmetric (LRS) cosmological models of Bianchi 
type-II representing clouds of geometrical as well as massive strings. Wang [T3] studied the Letelier model in the 
context of LRS Bianchi type-II space-time. Recently, Pradhan et al. [HI [TB] and Amirhashchi and Zainuddin 
[T7] obtained LRS Bianchi type II cosmological models with perfect fluid distribution of matter and string dust, 
respectively. Belinchon |18l 119] studied Bianchi type-II space-time in connection with massive cosmic string and 
perfect fluid models with time varying constants under the self-similarity approach respectively. Recently, Tyagi 
and Sharma |20j have investigated string cosmological models in Bianchi type-II space-time. 

Motivated by the above discussions, in this paper, we have investigated a new class of Bianchi type-II 
cosmological models for a cloud of strings by using the law of variation for generalized mean Hubble's parameter. 
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This approach is different from what the other authors have adapted. The paper is organized as follows. The 
metric and the field equations are presented in Section 2. Section 3 deals with exact solutions of the field 
equations with cloud of strings. Physical behavior of the derived model is elaborated in detail. Finally, in 
Section 4, concluding remarks are given. 

2 The metric and field equations 

We consider totally anisotropic Bianchi type-II line element, given by 

ds 2 = -dt 2 + A 2 [dx - zdyf + B 2 dy 2 + C 2 dz 2 , (1) 

where the metric potentials A, B and C are functions of t alone. This ensures that the model is spatially 
homogeneous. 

The Einstein's field equations ( in gravitational units c = 1, 8irG = 1) read as 

G) = -T], (2) 

where Gj — Rj — \Rg) is the Einstein tensor. The energy-momentum tensor Tj for a cloud of massive strings 
and perfect fluid distribution is taken as 

Tj = (p + p)v l vj + pg* - \x l Xj , (3) 

where p is the isotropic pressure; p is the proper energy density for a cloud strings with particles attached to 
them; A is the string tension density; v l — (0, 0, 0, 1) is the four-velocity of the particles, and x l is a unit space-like 
vector representing the direction of string. The vectors v % and x % satisfy the conditions 

v l v l = —x l x l = -1, v l x t — 0. (4) 

Choosing x l parallel to d/dz, we have 

x i = (0,0,C -1 ,0). (5) 

Here the cosmic string has been directed along z-direction in order to satisfy the condition = T 2 . As a 
result, the off-diagonal component of Einstein tensor, viz., G\ — z(T 2 — T^) vanishes. A detailed analysis about 
the choice of energy- momentum tensor for Bianchi type-II models is given by Saha |21] . 

If the particle density of the configuration is denoted by p p , then 

p = p p + A. (6) 

The Einstein's field equations ([2]) for line element ([lj with energy- momentum tensor ((3]), lead to the following 
set of independent differential equations: 

n i>r< q a2 

- -P , (7) 

= ~P , (8) 
-P + A, (9) 
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Here, and in what follows, an over dot indicates ordinary differentiation with respect to t. The energy 
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conservation equation T u = 0, leads to the following expression: 



which is a consequence of the field equations (j7|)-(TT0 
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3 Solutions of the Field Equations 



Equations (|T1)- (|10[) are four equations in six unknown parameters A, B, C, p, p and A. Two additional constraints 
relating these parameters are required to obtain explicit solutions of the system. 

First, we utilize the special law of variation for the Hubble's parameter given by Berman [22] , which yields 
a constant value of deceleration parameter. Here, the law reads as 



H = £(ABC)- 



(12) 



where £ > and n > are constants. Such type of relations have already been considered by Berman and 
Gomide [23] for solving FRW models. Later on, many authors (see, Kumar and Singh [24], Akarsu and Kilinc 
[25] and references therein) have studied flat FRW and Bianchi type models by using the special law for Hubble's 
parameter that yields constant value of deceleration parameter. 

Considering (ABC) ' as the average scale factor of the anisotropic Bianchi-II space-time, the average Hubble's 
parameter may be written as 



TJ I A B C 

h = tAa + b + c 



(13) 



Equating the right hand sides of ([12]) and (| 13[) . and integrating, we obtain 

ABC = (nit + ci)», (n^O) (14) 

ABC = c\e u \ (n = 0) (15) 

where C\ and c-i are constants of integration. Thus, the law (|12[) provides power-law (|14l) and exponential-law 
(fTS"]) of expansion of the Universe. 

Following Pradhan and Chouhan |26j . we assume that the component a\ of the shear tensor a" '. is propor- 
tional to the expansion scalar (8). This condition leads to the following relation between the metric potentials: 



A = (BC) m , 

where m is a positive constant. 

Now, subtracting ((5|) from ([7]) , and taking integral of the resulting equation two times, we get 



(16) 
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where C3 and C4 are constants of integration. In the following subsections, we discuss the string cosmology 
using the power-law (]14[) and exponential-law (]15j) of expansion of the Universe. 

3.1 String Cosmology with Power-law 

Solving the equations ([Tl]). (ITB]) and (flTj) . we obtain the metric functions as 



B(t) = a(n£t + ci)"< m + 1 ) exp 



1 3(l-m) 

C(t) = c A \ntt + ci) exp 



A(t) = (nlt + ci)< m +V 
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where M = (9m — 3 + mn + n)(3m — 3 + ran + n) and 
In the special case n = 3, we have 

= (3€t + ci)^r 



3^^+03(^+1) 
B(t) = a(3£t + ci) 3 '(™+ 1 ) exp 



(m + l) s 
144^ 2 m 2 
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Thus, the metric (JXJ) is completely determined. 

The expressions for the isotropic pressure (p), the proper energy density (p), the string tension (A) and the 
particle density (p p ) for the above model are obtained as 

3£ 2 [n(m + 1) - 3(m 2 - m + I)} . flj ._ 2 3(m + l)(n+ 1) , „ ± ad-") 

v = — — t ^ -(n£t + ci) 2 + '—(n£t + a) 

y (m + 1) 2 V ly 4(9m- 3 + mn + n) K ' 

3c 3 £(l-2m), ^ .a i (m + 1) 2 , . r iy-m) 2 i 

to + 1 r(9m — o + mn + rt) z 

- ,, q 2C3 ^ +1) + + cO-t™-] - cg(nft + *)-«, (24) 

t!(9to — 3 + mn + n) 



9t? 2 to(2 — to) . „ . 9 15 — 33771 — mn — n , „ , e(i-m) 
7 \_ n2 K nft + ci " +175 — ^(nft + ci "<- +1 > 

(TO + l) z 4(9777 — O + 77777, + 77) 

3C3£(1-2TO). , ,31 (TO+1) 2 , „ r l2(l-m) 1 

- ,, q 2C3 t +1 \ A rUt + c,)-^- 1 ! ^ c 2 (t^ + Cl )-, (25) 
£(9to — 3 + mn + n) 

A = M2i2m m lf ~ n) {nit + c l} - 2 + 2(nft + cQ-^ , (26) 



3£ 2 [3to(2-to)) + (3-?i)(1-2to)(to+1)1 / „ , 9 3(13 - 35777 - 377777 - 377), „ , 6(l-m) 

1 ^ / — ^(ntt + Cl ) - 2 + ^ I75 q— T- -^(rdt + c) 

(m + l) z 4(9m — 6 + mn + n) 



3c 3 l(l - 2m) 3_ x (m + 1) 2 [itosi_ 2 l 

+ TO+1 (T ^ + Cl) " " ^(9TO-3 + TO77 + 77) 2 ( ^ + Cl) ^ ' + > J 

- m 2C f: +1) + , (^ + c 1 )-[^- 1 ]-c 2 (77^ + Cl )^. (27) 

£(9?77 — O + TO77 + 77) 

The above solutions satisfy the energy conservation equation (TTTj) identically, as expected. 




Figure 1: Scale factors vs time with n = 0.5, to = 0.4, £ = 2, c\ = c 3 = 1 and C4 = 2. 

We observe that all the parameters diverge at t = —c\jn£. Therefore, the model has a singularity at 
t = —ci/n£, which can be shifted to t — by choosing c\ = 0. This singularity is of Point Type as all the 
scale factors vanish at t — —c\/n£. The cosmological evolution of Bianchi-II space-time is expansionary since 
all the scale factors monotonically increase with time (see, Fig.l). So, the Universe starts expanding with a big 
bang singularity in the derived model. The parameters p, p, p p and A start off with extremely large values. In 
particular, the large values of p p and A in the beginning suggest that strings dominate the early Universe. For 
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sufficiently large times, p p and A become negligible. Therefore, the strings disappear from the Universe for larger 
times. That is why, the strings are not observable in the present Universe. 
The rates of expansion in the direction of x, y and z are given by 

h x = - = -—-(net + a r 1 , 28 

A to + 1 

ff « - ! = ^r<- a + «>-' + e iSm -T + L + ,) ^ + + +«>-*. <»> 

H. 4 = ^+^<"« + <=■)-- g7q T M l, + e,)-SSi« - c s („« + «)-*. (30) 

O to+1 £(9to — 6 + mn + n) 

The average Hubble's parameter, expansion scalar and shear of the model are, respectively given by 

H=l(H x +H y + H z )=t{n£t + c 1 )- 1 , (31) 
6 = 3H = 3£{n£t + Ci)- 1 , (32) 

<J 2 = \[{H x -H y ) 2 + {H y -H z f + {H z -H x ) 2 ] 

3£ 2 (2m-l) 2 „ , 2 15-33TO-TOn-n, „ , 

= r xn2 nft + c i) - 775 =TI — ^( ntt + Cl ) 

(m + l) z 4(9to — 3 + mn + n) 

3c 3 £(l-2m). ^ .3 , (to + 1) 2 , . , \ mi- m) i 

to + 1 z z (9m — 3 + mn + n) z 



+ „ Q 2C3 f + 1) , (nft + cO-t™- 1 ! + c 2 Kt + 
£(9m-3 + mn + n) v dv y 



(33) 



The spatial volume (V) and anisotropy parameter (A) are found to be 

V = ABC = {ntt + c x )n, (34) 



7 2a 2 2(2m-l) 2 15 - 33m- tow -n , „ 6d-™) , 2 

A = = — (ntt + cA "0+i) + 

3iJ 2 (to+1) 2 6^ 2 (9TO-3 + TOn + n) v ; 

2c 3 (l-2m),_ Al , _ 3 +1 , 2(to + 1) 2 , _ 



*(m+l) + +3^(9TO-3 + TOn + n) 2( ^ + Cir 

4c3(m + 1} [nit + Cl )-t^- 3 ] + 2 An£t + Cl )- + 2 . (35 ) 



3£ 3 (9m-3 + mn + n) v ' U 2 

The value of DP (q) is found to be 



-l + -2=n-l, (36) 



H 
IT 2 

which is a constant. A positive sign of q, i.e., n > 1 corresponds to the standard decelerating model whereas 
the negative sign of q, i.e., < n < 1 indicates acceleration. The expansion of the Universe at a constant rate 
corresponds to n = 1, i.e., g = 0. Also, recent observations of SN la [37]- [33] reveal that the present Universe is 
accelerating and value of DP lies somewhere in the range — 1 < q < 0. It follows that in the derived model, one 
can choose the values of DP consistent with the observations. 

From the above results, it can be seen that the spatial volume is zero at t = —Cx/n£, and it increases with 
the cosmic time t. The parameters H x , H y , H z , H, 9 and a diverge at the initial singularity. These parameters 
decrease with the evolution of Universe, and finally drop to zero at late times provided m < 1. The mean 
anisotropy parameter asymptotically approaches to 2 ^^~^1 for n < 3 and to < 1. Thus, the dynamics of the 
mean anisotropy parameter depends on the values of n and m. The model does not approach isotropy provided 
to =/= 0.5 as may be observed from Fig. 1. In case to = 0.5, at late times, the directional scale factors vary as 

A(t) w (n£t + cx)n, B(t) » (ntt + a)*, C{t) » {nit + ci)» . 
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Therefore, isotropy is achieved in the derived model for m — 0.5. Since the present-day Universe is isotropic, 
we consider m = 0.5 in the remaining discussion of the model. 

For n < 1, the model is accelerating whereas for n > 1 it goes to decelerating phase. In what follows, we 
compare the two modes of evolution through graphical analysis of various parameters. We have chosen n = 2, 
i.e., q = 1 to describe the decelerating phase while the accelerating mode has been accounted by choosing n = 0.5, 
i.e., q — —0.5. The other constants are chosen as I = 6, ci = 1, c 3 — 1, m — 0.5. 




Figure 2: Pressure vs time. 
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Figure 3: Rest energy density vs time. 

Fig. 2 depicts the variation of pressure versus time in the two modes of evolution of the Universe. We observe 
that the pressure is positive in the decelerating Universe which decreases with the evolution of the Universe. But 
in the accelerating phase, negative pressure dominates the Universe, as expected. In both cases, the pressure 
becomes negligible at late times. 

The rest energy density has been graphed versus time in Fig. 3. It is evident that the rest energy density 
remains positive in both modes of evolution. However, it decreases more sharply with the cosmic time in the 
decelerating Universe. 

Fig. 4 and Fig. 5 show the behavior of particle energy density and string tension versus time in the 
decelerating and accelerating modes, respectively. We see that p p > A, i.e., the particle energy density remains 
larger than the string tension density during the cosmic expansion , especially in early Universe. This shows 
that massive strings dominate the early Universe (see, Refs. [UQjj]). Further, it is observed that for sufficiently 
large times, p p and A tend to zero. Therefore, the strings disappear from the Universe at late times. 

According to Ref. |19j . since there is no direct evidence of strings in the present-day Universe, we are in 
general, interested in constructing models of a Universe that evolves purely from the era dominated by either 
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Figure 4: Particle energy density and string tension vs time for q = 1. 
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Figure 5: Particle energy density and string tension vs time for q = —0.5. 



geometric string or massive strings and ends up in a particle dominated era with or without remnants of strings. 
Therefore, the above model describes the evolution of the Universe consistent with the present-day observations. 
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Figure 6: p + p and p — p vs time for q = —0.5. 
From Fig. 3 to Fig. 6, we observe the following: 

00 9 > 

(ii) Pp > 

(iii) p + p > 

(iv) p — p > 0. 

This shows that the weak and dominant energy conditions are satisfied in the derived model. 

3.2 String Cosmology with Exponential-law 

Solving the equations ([IS]), (fTS]) and (fTT|). we obtain the metric functions as 

3m£ 



A(t) 



exp 



m + 1 



B(t) — CiC^ +1 exp 
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c, exp 
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(37) 
(38) 
(39) 



The expressions for the isotropic pressure, the proper energy density, the string tension and the particle 
density for the derived model are obtained as 



P = 



9£ (m 2 ~m + l) c 2 m+1 (m + 1) 6i(m-i) f 2c^c 2 m+1 (m + 1) mrn-a) 
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6(m-l) 3(m-3) 

9£ 2 m(2 - m) c 2 m+1 (5 - 11m) w(m-o t 2c 3 c 2 m+1 (m + 1) s<(m-3) f 
9 = (m + 1) 2 + 4(3m - 1) 6 ™ +1 3^(3m - 1) 6 

12(m-l) 

3c 3 (l-2m) m _ c 2 m+1 (m + 1) 2 12^(^-1) t _ c\ m 
c|(m+l) 9^ 2 (3m-l) 2 c| 1 j 

x 9i? 2 (2m-l) Hfeii 6<(m-i) . 

A = 1 — ' + 2c 2 m+1 e^+^', 42 

m+1 

6(m-l) 3(m-3) 

9£ 2 (l + m-3m 2 ) c 2 m+1 (13 - 25m) s<(m-i) f 2c 3 c 2 m+1 (m + 1) 3^-3) t 
Pp (m + 1) 2 + 4(3m- 1) 6 m+ 3^(3m - 1) 6 

12(m-l) 



3c3(l-2m) _ c 2 m+1 (m+1 ) 2 u t c§ ^_ 6tt 
cf(m+l) 9^ 2 (3m-l) 2 



The energy conservation equation (fTTj) is satisfied identically by the above solutions, as expected. 
The directional Hubblc's parameters are given by 

H. = -^r, (44) 
m + 1 

6(m-l) 



_ -~ ^2 (m+i) «fe^a t C3 3^ 
^ _ m+l + 3£(3m-l) 6 + cf ' (45) 

6(m-l) 

_ 3f(l-m) c 2 "' +1 (m + 1) M^l> t c 3 3ft 
2 " m+1 3*(3ro-l) 6 cf ' 1 ° j 

Hence the average generalized Hubble's parameter is given by 

H = i. (47) 

From equations (|44l) - (|4"Tj) . we observe that the directional Hubble's parameters are time dependent while the 
average Hubble's parameter is constant. 

The expressions for kinematical parameters, i.e., the scalar of expansion, shear scalar, the spatial volume, 
average anisotropy parameter and deceleration parameter for the derived model are given by 

6 = 31, (48) 



6(m-l) 3(m-3) 

2 _ 3^ 2 (2m-l) 2 c 2 m+1 (5 -11m) <u(m-i) f 2c 3 c 2 m+1 (m + 1) 3<(m- 3) , 
(m + 1) 2 4(3m- 1) 6 ™ +1 + 3€(3m - 1) 6 

3c3(l-2m) 3^ c 2 - +1 (m + 1) 2 ™im=H t c\ m 
4(m+l) + 9£ 2 (3m-l) 2 c§ ' K ' 

V = c\e u \ (50) 

6(m-l) 3(m-3) 



j 2(2m-l) 2 c 2 m+1 (5 -11m) 6*(m-p f 4c 3 c 2 m+1 (m + 1) aijm-s) t 

(m + 1) 2 6£ 2 (3m - 1) 6 ™ +1 + 9F(3m- 1) 6 

12(m-l) 

2c 3 (l — 2m) _ v/ 2c 2 m+1 (m+1) 2 i2i(m-i) . 2c 2 _ fi „ 

-e + ;L M , rry^e ™+i +^%e 6W , (51) 



£ 2 cf(m+l) 27^ 4 (3m-l) 2 3£ 2 



2 



g=-l. (52) 

Recent observations of SN la [37]- [33] suggest that the Universe is accelerating in its present state of evolution. 
It is believed that the way Universe is accelerating presently; it will expand at the fastest possible rate in future 
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and forever. For n = 0, we get q = — 1 ; incidentally this value of DP leads to dH/dt = 0, which implies the 
greatest value of Hubble's parameter and the fastest rate of expansion of the Universe. Therefore, the derived 
model can be utilized to describe the dynamics of the late time evolution of the actual Universe. So, in what 
follows, we emphasize upon the late time behavior of the derived model. At late times, we find 



In particular, for to — \, we have 



9e 2 (m 2 -m + 1) 
(to + l) 2 

^ 9£ 2 m(2 - to) 
" (m + 1) 2 ' 

^ 9£ 2 (2m- 1) 

TO + 1 

9^ 2 (1 + to-3to 2 ) 
(m + 1) 2 

r2 _ 3 l 2 (2m- l ) 2 
(to + l) 2 

2(2m- l) 2 



(53) 
(54) 
(55) 



2 ^ - -l; 



(m + 1) 2 



P = -P, 
A w 0, 
P P « 3^ 2 , 
a 2 «0, 



(58) 



This shows that vacuum energy dominates the Universe at late times, which is consistent with the ob- 
servations. Strings disappear and the Universe evolves with constant particle energy density. The shear and 
anisotropy parameter become negligible. So the Universe becomes isotropic. 



4 Concluding Remarks 

In this paper, a spatially homogeneous and anisotropic Bianchi-II space-time representing massive strings in 
general relativity has been studied. The main features of the work are as follows: 

• The models are based on exact solutions of the Einstein's field equations for the anisotropic Bianchi-II 
space-time filled with massive strings. 

• The singular model (n ^ 0) seems to describe the dynamics of Universe from big bang to the present epoch 
while the non-singular model (n = 0) seems reasonable to project dynamics of future Universe. 

• In the derived models, m=\ turns out to be a condition of isotropy. 

• In the present models, the weak and dominant energy conditions are satisfied, which in turn imply that 
the derived models are physically realistic. 

• The singular model presents the dynamics of strings in the accelerating and decelerating modes of evolution 
of the Universe. It has been found that massive strings dominate the Universe, which eventually disappear 
from the Universe for sufficiently large times. This is in agreement with the astronomical observations. 

• The non-singular model for m = \ predicts a Universe dominated by vacuum energy, which is consistent 
with the predictions of current observations. 
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